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We study concentrated binary colloidal suspensions, a model system which has a glass transition as the volume fraction (j) of 
particles is increased. We use confocal microscopy to directly observe particle motion within dense samples with (]) ranging 
from 0.4 to 0.7. Our binary mixtures have a particle diameter ratio d^jd^ = 1/1.3 and particle number ratio N$/Nl = 1.56, 
which are chosen to inhibit crystallization and enable long-time observations. Near the glass transition we find that particle 
dynamics are heterogeneous in both space and time. The most mobile particles occur in spatially localized groups. The length 
scales characterizing these mobile regions grow slightly as the glass transition is approached, with the largest length scales seen 
being ~ 4 small particle diameters. We also study temporal fluctuations using the dynamic susceptibility %4, and find that the 
fluctuations grow as the glass transition is approached. Analysis of both spatial and temporal dynamical heterogeneity show that 
the smaller species play an important role in facilitating particle rearrangements. The glass transition in our sample occurs at 
§ g « 0.58, with characteristic signs of aging observed for all samples with (]) > § g . 
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.1 Introduction 

As the temperature of a glass-forming liquid is lowered, the 
viscosity rises by many orders of magnitude, becoming exper- 
imentally difficult to measure, with little change in the struc- 
ture^ - —. The origin of the slowing dynamics is not yet clear, 
despite much prior work. One intriguing observation is that 
■as a sample approaches the glass transition, the motion within 
the sample becomes spatially heterogeneous 4 ' 5 . While overall 
motion within the sample slows, some regions exhibit faster 
dynamics than the rest, and over time these mobile regions ap- 
pear and disappear throughout the sample. Particles within the 
mobile region move cooperatively, forming spatially extended 
clusters and strings 6 . 

One technique for studying the glass transition is the use of 
colloidal suspensions 7 . These are composed of small solid 
particles suspended in a solvent. The particles need to be 
small enough to undergo Brownian motion, so particle diam- 
eters are typically 10 — 5000 nm. The key control param- 
eter is the volume fraction (j). For a monodisperse sample 
(all particles similar in size), the sample becomes glassy for 
<|) >ty g « 0.58-^. The glass transition in colloidal samples has 
been studied extensively by light scattering, microscopy, and 
other techniques. Colloidal samples exhibit many behaviors 
seen in molecular glasses, such as dramatic increases in vis- 
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croscopic disorder 17 , spatially heterogeneous dynamics-^—, 
aging behavior for glassy samplesi^22r— , and sensitivity to 
finite size effects 28 - 29 . Light scattering allows careful study of 
the average behavior of millions of colloidal particles, while 
microscopy techniques observe the detailed behavior of a few 
thousand particles. These complementary techniques have re- 
sulted in connections between different aspects of glassy be- 
havior: for example, showing that aging is temporally and 
spatially heterogeneou s 22 i 23 , and connecting dynamical het- 
erogeneity with the slowing relaxation time scale s 8 i 13 ' 19 i 20 i 30 . 

In this paper, we study the glass transition of binary col- 
loidal suspensions using confocal microscopy. We use binary 
suspensions (mixtures of two particle sizes) to inhibit crystal- 
lization. This allows us to take data over many hours, a time 
scale in which a monodisperse sample would crystalliz e 31 i 32 . 
Furthermore, this lets us investigate the role the two particle 
species play in the dynamics; prior work has suggested that 
small particles play a lubricating role in the local dynamics 33 . 
Prior studies have also seen a connection between the local 
structure and the mobility of particles^"—. Using a binary 
sample results in more obvious structural variations due to 
spatial variability of the composition, helping highlight how 
structure influences the dynamics. 

The confocal microscope enables direct visualization of the 
interior of the sample, and we follow the motion of several 
thousand colloidal particles within each sample 38 . Particles 
move in spatially heterogeneous groups, and we character- 
ize this motion using two-particle two-time correlation func- 
that have previously been used on monodisperse 
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suspensions^. From these we extract a length scale for the 
heterogeneity, which increases as the glass transition is ap- 
proached. By simultaneously tracking both large and small 
particles, we can observe the similarities and differences be- 
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tween the two species' dynamics. In particular, we see that 
small local composition fluctuations influence mobility. As 
might be expected, regions with more large particles tend to be 
less mobile, while those with more small particles tend to be 
more mobile. Additionally, we study temporal heterogeneity 
using a different correlation function, the dynamic susceptibil- 
ity 2£ 4 i2"d£, which has not been previously applied to colloidal 
data. As measured by this correlation function, the temporal 
heterogeneity increases as the glass transition is approached. 

2 Experimental Method 

We prepare suspensions of poly-(methyl-methacrylate) 
(PMMA) colloids stabilized sterically by a thin layer of 
poly-12-hydroxystearic acid 7 . We use a binary mixture with 
a large particle mean radius cil = 1 .55 /jm and small particle 
mean radius as = 1.18 /jm, the same size particles used in a 
prior study by our group 28 . The polydispersity is 5%; each 
individual particle species can crystallize in a monodisperse 
suspension. Separately from the polydispersity, the mean par- 
ticle radii each have an uncertainty of ±0.02 /jm. The number 
ratio of small particles to large particles is Ns/Nl = 1.56, 
resulting in a volume fraction ratio 0s/0£ ~ 0.70. The control 
parameter is the total volume fraction = 0$ + 0^. Crystal- 
lization and segregation were not observed to occur during the 
course of our measurements. All particles are fluorescently 
dyed and suspended in a density- and index-matched mixture 
of decalin and cyclohexyl bromide to prevent sedimentation 
and allow us to see into the sample. Particles are slightly 
charged as a result of the dyeing process and this particular 
solvent mixture 32 . Nonetheless, we use the hard sphere 
volume fraction as the control parameter. The hard sphere 
radii (ai,as) are determined from diffusion measurements of 
the individual species at dilute concentrations (0 < 0.01). 

Suspensions are sealed in microscope chambers and con- 
focal microscopy is used to observe the particle dynamics at 
ambient temperature 38 ' 46 . A representative two-dimensional 
image is shown in Fig.Q] A volume of 55 x 55 x 20 /jm 3 can 
be taken at speeds of up to 1 Hz. (As will be shown later, 
in these concentrated samples, particles do not move signifi- 
cantly on this time scale.) To avoid influences from the walls, 
we focus at least 25 /jm away from the coverslip. 

Within each three-dimensional image, we identify both 
large and small particles. This is accomplished with a single 
convolution that identifies spherical, bright regions 47 ; the con- 
volution kernel is a three-dimensional Gaussian with a width 
chosen to match the size of the image of a large particle. Each 
local maximum after the convolution is identified as a parti- 
cle^ 7 -. The distribution of particle brightnesses is bimodal with 
little overlap, and so small and large particles can be easily dis- 
tinguished. Our method is the same as is often used to measure 
particle positions in two dimensions, which normally achieves 
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Fig. 1 A two-dimensional image of our sample taken by a confocal 
microscope. The scale bar represents 10 /jm. 

sub-pixel resolution in particle positions 47 . However, given 
that a single convolution kernel is used to identify both particle 
types, when applied to our binary samples we do not achieve 
sub-pixel accuracy. Instead, our uncertainty in locating parti- 
cle positions is linked to the pixel size and is 0.2 /jm in x and y, 
and 0.3 /jm in z- We do have accurate discrimination between 
large and small particles with this method, with less than 1 % 
of the particles misidentified, checked by visual inspection. 
For a few particles, it is hard to distinguish if they are large 
or small (because they are small but unusually bright, or large 
but unusually dim). These particles are assigned to the size 
that they appear to be the majority of the time. 

After identifying the particle positions, they are tracked us- 
ing standard softwar e 46 ' 47 . The key requirement is that parti- 
cles move less between time steps than their interparticle spac- 
ing, which is easily satisfied in our dense glassy samples. We 
take images once every 10 - 150 s, depending on the volume 
fraction, in each case making sure that the acquisition rate is 
sufficiently rapid to capture all particle movements. 

3 Results and Discussion 
3.1 Structural characteristics 

We begin by looking at the structure of the binary sample. 
Shown in Fig.|2]is the pair correlation function g(r) of a sam- 
ple with volume fraction = 0.57. g(r) relates to the likeli- 
hood of finding a particle a distance r away from a reference 
particle. The three curves correspond to small-small, small- 



large, and large-large particle pairs, and their first peak posi- 
tions are at approximately 2as = 2.36 /jm, as + «l = 2.73 /jm, 
and 2ai = 3.10 pm. The first peaks of the three curves are 
fairly broad due to the uncertainty in locating particle posi- 
tions, and also in part due to the particle polydispersity. 
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Fig. 2 The pair correlation function g(r) for a sample with volume 
fraction <|> = 0.57. The different curves are for small-small pairs, 
small-large pairs, and large-large pairs, as indicated. Due to particle 
tracking uncertainties (±0.2 /jm in x and y, ±0.3 /jm in z), the 
measured separation r between any pairs of particles has an 
uncertainty of ±0.5 jum, which significantly broadens the peaks of 
g(r) and diminishes their height. 



3.2 Dynamical slowing 

We wish to show how the motion of particles slows as the 
volume fraction increases and approaches the glass transi- 
tion. Figure [3] shows results of the mean square displacement 
(MSD) (Arf) of large and small particles, where A?,- = Ar,(Af) 
denotes the displacement of 2-th particle in lag time At, and 
the brackets indicate an average over all particles and times 
observed. Figure [3] shows that as the volume fraction in- 
creases, particle motion slows significantly, as expected. At 
(|) = 0.4, small particles take tens of seconds to move a dis- 
tance a c = 1.4 jum 2 ; at (j) = 0.54 the time has grown to more 
than 10 4 s. For the lowest volume fraction samples, comparing 
the two particle species, we find that (Ar?) / (Arj) w a^/as, as 
expected from the Stokes-Einstein-Sutherland equation 48 ! 49 . 

As the volume fraction increases, the MSD plots show a 
characteristic "cage trapping" plateau. Particles cannot dif- 
fuse freely, but instead are "caged" by their nearest neigh- 
bor particle s 34 i 5Q ~ ^l. The upturn in the MSD curve is iden- 
tified with rearrangements of the cage, allowing the particle 
to move to a new location, perhaps caged by different par- 
ticles. Although the smaller particles diffuse faster than the 



large particles, MSD curves for both show upturns at similar 
time scales, indicating that their dynamics are strongly cou- 
pled—. Note that the height of the plateaus (0.1 - 0.2 pm 2 ) 
are larger than in prior work~ , because of the larger particle 
tracking noise present in this binary experiment as compared 
to the prior work with monodisperse particles. The noise re- 
sults in apparent displacements (Ax 2 ) = (Ay 2 ) — (0.2 /jm) 2 , in- 
dependent of At, and so the minimum reliable value for Fig. [3] 
is limited to ~ 0.08 /jm 2 . In particular the difference between 
the (j) = 0.58 and (j) = 0.66 data is probably not significant, but 
due to higher tracking errors for the higher volume fraction 
data. The noise prohibits careful analysis of the mean square 
displacement data along the lines of prior wor k 55 ' 56 . 
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Fig. 3 (Color online) A log-log plot of mean square displacement 
versus time lag for large particles (a) and small particles (b). Note 
that our resolution means we cannot accurately measure mean 
square displacement values less than 0.1 /jm 2 , and thus the plateau 
height of curves for the highest volume fraction data is set by this 
limit, rather than the dynamics. However, the slight upturn for those 
curves at large values of At is above our resolution limit and thus 
real. Due to noise in (At}), the data shown are (Ax 2 + Ay 2 }. 

For the samples with (]) > (j) ? w 0.58, the MSD curves are 
nearly flat, suggesting that on our experimental time scales, 
these samples behave as glasses. Glasses are fundamen- 
tally non-equilibrium systems, so that physical properties for 
glasses depend on the preparation history in general and, in 
particular, the time since they were initially formed. This 
time-dependence is known as aging, and can be quantified by 
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Fig. 4 (Color online) MSD plot for the large particles from a 
(|> = 0.59 sample. The three curves correspond to time regimes with 
< t < 4800 s, 4800 < t < 9600 s, and 9600 < t < 14400 s. Note 
that in the first time regime, the sample ages appreciably, so the 
"early" mean square displacement data should be interpreted with 
caution. The time-dependence of (Ar 2 ) is clearly seen, indicating 
the presence of aging. As = 0.59 is the lowest volume fraction in 
which this behavior is seen, we conclude that the glass transition 
occurs at § g ~ 0.58. Due to noise in (Az 2 ), the data shown are 
(Ar 2 ) = (Ax 2 + Ay 2 ). 



examining the MSD at different times since the start of the 
experiment 22 i 26,33 . Figure [4] shows MSD data from (j) = 0.59. 
The trajectory data are broken into three equal duration seg- 
ments and the MSD calculated within each segment; for the 
older segments, the MSD curve decreases in height (see the 
caption for details). The sample is most active immediately af- 
ter being formed, and continues to slow down as time elapses. 
The aging of the MSD appears in samples for > 0.59, while 
no samples for (j) < 0.58 show aging. From the onset of aging, 
we conclude that the glass transition point is at volume frac- 
tion (j) w 0.59, similar to that seen for monodisperse samples 7 . 
Note that our particle size uncertainty of ±0.02 /jm (radius) 
leads to a systematic volume fraction uncertainty, so our esti- 
mate is (j) g = 0.59 ± 0.02 as a comparison with other work. 

Particles involved in a cage rearrangement event move sig- 
nificant distances compared to when they are caged, and prior 
work noted that the distribution of displacements is unusu- 
ally broad on the time scale of the rearrangement 2aiZ. This is 
quantified by calculating the non-Gaussian parameter ot2(Af ), 
which for a one-dimensional distribution of displacements is 
defined as 

(Ax 4 ) 



a 2 (At) 



1. 
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3 (Ax 2 ) 

where Ax = x(f + At) —x(t) denotes the x displacement for 
time lag Af 58 , and the angle brackets indicate an average over 
all particles and all initial times t. If the distribution of dis- 
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Fig. 5 (Color online) A semi-log plot of the non-Gaussian parameter 
0C2 versus lag time for large particles (a) and small particles (b). The 
curves for = 0.40,0.42 are essentially indistinguishable. 



placements Ax is Gaussian, then 0C2 = by construction. If 
events with large displacements are more common than would 
be expected from a Gaussian distribution, then (X2 > 0. Fig- 
ure [5] shows the results of the non-Gaussian parameter (NGP) 
for one-dimension displacements of large and small particles. 
The curves peak at time scales where cage rearrangements are 
most important 6 - 20 , and thus coincides with the upturn of the 
MSD curves. 

(X2 is fairly sensitive to experimental noise, although fortu- 
nately the "signal" this parameter measures comes from parti- 
cles moving large distances, which are less susceptible to par- 
ticle position uncertainty. The differences between the data 
shown in Fig. |5]and the data of Ref. 20 are probably more due 
to noise and uncertainties in measuring (j), rather than system- 
atic differences between our binary experiments and the prior 
monodisperse experiments 20 . The peak heights of (X2 for our 
small particles are similar to those seen in Ref M, for roughly 
similar volume fractions. One notable difference is that Ref. 20 
found that for glassy samples, 0C2 started high (1-2) and de- 
creased steadily, whereas in our current data, 0C2 remains fairly 
small for all At, with unclear dependence on At. This is proba- 
bly due to our noise, as for glassy samples, few particles move 
distances large enough to be outside of our particle tracking 
uncertainty. 
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Figure [5] also reveals that the motions of the small particles 
are more dynamically heterogeneous, with the maximum NGP 
peaking above 1.5 for the small species but only reaching 0.8 
for the large species. This is consistent with recent observa- 
tions of aging binary colloidal glasses, which likewise found 
the small particles had more non-Gaussian motion—. 

From Figs. [3] and [5] we conclude that the dynamics of the 
large and small particles are qualitatively the same, although 
with small quantitative differences. In particular, the time 
scale over which particles escape cages is the same for both, 
as is the time of peak non-Gaussianity. In much of the subse- 
quent analysis, therefore, we consider both species together in 
order to obtain better statistical validity. 



3.3 Local environment influences mobility 

We wish to understand the origins of dynamical heterogene- 
ity. For a hard-sphere system, or an overdamped system such 
as our experimental colloidal suspension, the only variable is 
the local structure. Clearly structure has some relation with 
particle mobility 36 - 37 , although this relationship may be diffi- 
cult to see and not directly predictive in nature 35 . Prior work 
found that more disordered environments are weakly corre- 
lated with higher particle mobility 26 ' 34 , and a recent study of 
aging binary colloidal glasses found a relation between the lo- 
cal composition and the mobility—. 

We quantify a particle's local environment by counting its 
nearest neighbors Nnn, defined as particles closer than the first 
minimum of the pair correlation function for the large parti- 
cles, 4.1 /jm (Fig.|2]i, and distinguish between large and small 
neighbors. Figure [6] shows that the number of neighbors of 
a given type has a strong influence on the mobility of a par- 
ticle. Particles with more large neighbors have, on average, a 
lower mobility, while those with more small neighbors a larger 
mobility. These observations agree with studies of aging in bi- 
nary colloidal glasses, 33 and are reminiscent of prior Theolog- 
ical observations of binary suspensions 59-61 which noted that 
binary mixtures have lower viscosities than single-component 
samples with equivalent total volume fraction. The reason- 
ing is that binary suspensions can in general be packed to 
higher volume fractions, and so have more free volume than 
monodisperse samples at the same volume fraction. Figure [6] 
suggests that the small particles indeed "lubricate" large par- 
ticles, as previously proposed 61 . Conversely, large neighbors 
significantly inhibit the motion of both large and small parti- 
cles. The lubricant effect for large particles (which have less 
free volume) is less pronounced, agreeing with prior observa- 
tions of monodisperse suspensions which found that regions 
of mobility correlated with regions of larger free volume 34 . 
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Fig. 6 Large particle mobility as a function of the number of large 
and small nearest neighbors Nnn- The panel (b) shows small 
particle mobility. Mobility is very sensitive to the number of large 
neighbors, decreasing sharply as the number of large neighbors 
increases. The number of smaller neighbors has a weaker, yet 
measurable, impact. These data are for volume fraction (]) = 0.53, 
using a time scale At = 3780 s to define displacements. 



3.4 Cooperative motions 

Prior work has shown that the higher mobility molecules in a 
supercooled liquid are distributed in a spatially heterogeneous 
fashion 5,6 . In monodisperse colloidal systems, direct imaging 
using microscopy found that particles rearrange in coopera- 
tive groups 18 " 20 . Following the prior work, we characterize 
the cooperative nature of colloidal rearrangements by study- 
ing the dynamics over a time scale At* that corresponds to the 
maximum of the NGP^S. The maximum displacement of a 
particle over that time D, is defined as 



Di{t) := max (|r,-(f 2 ) -n{h)\) 
t,t+M* 



(2) 



where max rr+ Af* {X) is the maximum value of X using times 
t\, t2 such that t < t\ < ?2 < t + At*. Taking the maximum dis- 
placement results in a quantity that is less sensitive to random 
Brownian motion than the ordinary displacement Ar. Follow- 
ing prior work, 6 ' 20 a threshold D* ((])) is chosen such that on 
average, 5% of the particles at any given time have D,-(f ) > D* . 
These particles are termed "mobile particles" and generally 
are the ones undergoing cage rearrangements. (Note that at 
any particular time, the fraction of particles matching this def- 
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the vector and scalar spatial-temporal correlation functions 41 
S vec (R,At) and S sc i{R,At) defined as 
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Fig. 7 (Color online) Snapshots of system for = 0.54 (upper) and 
(]) = 0.66 (bottom). The red (dark gray) spheres are small mobile 
particles, and the blue (light gray) spheres are large mobile particles. 
Mobile particles are defined as those making the largest 
displacements at this particular moment in time; see text for further 
details. We set the time lag for the displacement as the cage 
breaking time scale (the peak time of the NGP ?ngp) which is 
At* = 3000 s for <b = 0.54 and t* = 4000 s for d> = 0.66. 



inition is not required to be 5%—.) 

Figure [7] shows snapshots of our system, highlighting the 
mobile particles. Clusters of these mobile particles are visible, 
in agreement with previous work which found similar mobile 
regions 6 ! 20 ' 33 . The clusters are somewhat smaller than those 
seen previously in single-component colloidal suspension 20 ; 
apparently the dynamics in binary mixtures are less spatially 
heterogeneous. Our result is in agreement with the results of a 
simulation study for polydisperse hard-disk systems 62 , which 
found that polydispersity reduces dynamic heterogeneity. It is 
also apparent in Fig. [7] that the small particles dominate the 
motions for the lower volume fraction (top, (j) — 0.54) whereas 
the two species contribute more equally in the glassy sample 
(bottom, (]) = 0.66). 



3.5 Length scales of spatial dynamical heterogeneity 

Pictures such as Fig. Q are qualitative evidence of dynami- 
cal heterogeneity. For quantitative information, we consider 
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The vector function S vec {R,At) characterizes correlations in 
the vector displacements Ar; = ry(r + At) — n(t); the sim- 
ilar function S sc i(R,At) uses the scalar displacement Sr, = 
| A7j | — ( | A7j | ) . The angle brackets () denote an average over all 
particle pairs with separation R at initial time t as well as an av- 
erage over t . The denominators of both correlation functions 
are averaged over all particles and time, and do not depend 
on R. The correlation function defined by eqn (01 indicates a 
vector correlation, and that defined by eqn (@]i a scalar correla- 
tion. If particles correlate perfectly, the correlation functions 
are unity. These correlation functions give information about 
spatial correlations for fixed lag time At, and about temporal 
dependence of the correlations for fixed separation R. We cal- 
culate these functions for all pairs of particles, without con- 
cern for the particle sizes, both to improve our statistics and 
because we do not find significant differences for large and 
small particles only. Most of the "signal" of correlated motion 
comes from the particles undergoing larger than average dis- 
placements, and so the results are less sensitive to the particle 
tracking uncertainties than the mean square displacement. 

Figure [8] shows the lag time dependence of these correla- 
tion functions, in which the distance R is set as the first peak 
distance of the small-large pair correlation function g(r) (the 
solid line in Fig. [2). At time scales larger than those shown in 
Fig. [8] the results become too uncertain, due to lack of data. 
In intermediate volume fraction region ((j) < (j) g ), both corre- 
lation functions increase with At. For the two lowest volume 
fractions ((j) = 0.42,0.49), the correlation functions eventually 
decrease at large Af , but our data do not extend to large enough 
At to see this for higher volume fractions. Overall, in con- 
junction with Fig. [3] Fig. [8] suggests that larger motions are 
more correlated with the motions of their neighboring parti- 
cles. This agrees with prior experiments 20,42 . The amplitude 
of the correlation decreases as (j) increases, with the exception 
of the (j) = 0.49,0.53 data which are similar. 

For glassy samples ((j) > 0.6) the correlation functions are 
small, suggesting that there is little correlation of the motion 
of neighbors. This is both because there is little overall mo- 
tion in glassy samples (see Fig. [3j and also the motion that 
does occur is dominated by Brownian motion within the cage, 
which is less correlated that the motions responsible for cage 
rearrangements 34 . Furthermore, it is probably erroneous to 
even consider time-averaged correlation functions for glassy 
samples, as in this current work the dynamics slow with time 
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Fig. 8 (Color online) Plot of the correlation functions in which the 
distance Rq = as + ai — 2.73 /jm is set at the sum of the particle 
radii, (a) represents the vector correlation [eqn ([3)], and (b) the 
scalar correlation [eqn ©]. 



(Fig. 2]i, and so a time average is of dubious validity. (Ear- 
lier work studied well-aged samples where the dynamics were 
only slow aging, and thus a time average was more sensi- 
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Fig. 9 Semi-log plot of the spatial correlation functions of = 0.53, 
where the time lag is set as A/ngp = 2000s. The solid line is 
5,. er (^,fNGp) and the dashed line is S sc i(R,fyfQp). The dotted line 
represents an exponential function with a decay length of 5.8 /jm, a 
good fit to both functions in this particular case. 
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Fig. 10 (Color online) The relationship between the length scales 
and the volume fraction. Shown are the length scales for the vector 
correlation function (open triangles) and the scalar correlation 
function (closed circles). The symbols indicate the average value, 
and the error bars show the range of values found for different lag 
times Af . These length scales are extracted from the correlation 
function for all particles (large + small). 



To consider the spatial dynamical heterogeneities, we plot 
the correlation functions as a function of R in Fig. [9] (for = 
0.54; results for other volume fractions are similar). For small 
separations around/? = 3.5 /jm, there is a dip in the correlation 
functions, which corresponds to the dip in the small-large pair 
correlation function at the same position (solid line in Fig.|2]i; 
the peak around R = 2.8 /jm likewise corresponds to the peak 
of the small-large pair correlation function. Thus, a particle's 
motion is correlated with that of its nearest neighbors, while 
particles separated by a less structurally favorable distance are 
less likely to have correlated motion. 

We fit our data with an exponential function S ~ 
Aexp(— R/Z,) and extract the decay length Figure [TOl shows 
both the vector (triangles) and scalar (circles) decay lengths as 
a function of the volume fraction. The length scales are essen- 
tially constant until close to ty g « 0.58, when they show a sharp 
increase. Our data are too noisy to draw conclusions about 
how the length scales grow near (j) ? , although simulations of 
binary Lennard-Jones mixtures did not find a divergence 63 . 
The largest length scale seen is rs 10 8as « 6.5az,, simi- 
lar to prior studies of monodisperse colloids 42 . For (j) > § g , the 
scalar length seems large and the vector length decreases, al- 
though as noted above, the data should be treated with caution 
as the samples are aging. 
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Fig. 11 (Color online) Surface plot of %4 for large particles within a 
sample with § = 0.52. 

3.6 Temporal dynamical heterogeneity 

The prior subsection showed that the motion of colloidal par- 
ticles in our dense samples are spatially heterogeneous. We 
now study their temporal heterogeneity, using the four-point 
susceptibility %4 which measures the correlation in dynam- 
ics between any two points in space within some time win- 
do w 43,45 . The actual value of %4 is a measure of the average 
number of particles whose dynamics are correlated, which in 
turn relates back to the spatial heterogeneity 14 . 

Here we only compute the self contribution to %4, since it 
has been shown to be the dominating ter m 43 ' 45 . The self part 
is computed from temporal fluctuations in particle mobility, 
where a particle is defined to be mobile if its displacement 
over some time interval At is larger than some threshold dis- 
tance AL^ 5 -. Using this definition, each particle at each time 
can be labeled mobile or immobile and the fraction of mobile 
particles Q(t) can be computed for each frame recorded. Q(t) 
varies from frame to frame due to the spatial heterogeneity. 
The temporal fluctuations in Q(t) are quantified by the self 
part to 5(4 and written as 

X4=N[(Q(tf) t -{Q(t))j], (5) 

where N is the number of particles. N also varies from frame 
to frame as particles move in and out of the field of view; we 
average N over all frames and use (N) in eqn (O. (The factor 
of N arises because the variance scales inversely with particle 
number.) Note that %4 measures temporal fluctuations in mo- 
bility without regard for the spatial correlations between mo- 
bile particles, whereas the correlation functions S vec and S sc i 
studied in the previous section measured spatial correlations 



of mobility without regard for the temporal correlations. In a 
sense, then, these two methods of analysis are complementary. 

From eqn (O it's evident that %4 will depend on our choice 
of AL and At as shown in Fig. QT| where %4 is plotted for the 
larger particles within a = 0.52 sample for various values 
of At and AL. This plot shows that 5(4 is characterized by a 
function that has a maximum at (Af max , AL max ). This maxi- 
mum in %4 indicates a typical timescale Af max where the dy- 
namics are most heterogeneous, and likewise AL max indicates 
a typical length scale distinguishing caged motions from cage 
rearrangements. 

Figure [T2l shows plots of %4(At,AL = AL max ) for the larger 
(a) and smaller (b) particles. The value of %4 is larger in mag- 
nitude for the smaller particles regardless of 0, demonstrating 
that the dynamics of the smaller particles are more temporally 
heterogeneous. Prior work by Lynch et al— showed a similar 
relative mobility; our results build upon this by showing that 
smaller particles also experience larger fluctuations, and thus 
exhibit more anomalous spatial and temporal behavior. We 
also see that %4 grows in amplitude as (j) increases, but then 
drops for the glassy sample (0 = 0.59), matching the results 
of the prior subsections where changes were seen at the glass 
transition. 

The plots in Fig.[T2]all show a maximum in %4 at a well de- 
fined Af max , and that Af max for the various volume fractions oc- 
cur at timescales close to where 0C2 shows a maximum in Fig. [5] 
and caging rearrangements become prominent (the "knee" in 
Fig. [3). The coincidence of maxima in 5(4 and 0C2 suggests 
that local cage rearrangements are the largest contributor to 
the temporal fluctuations. Since small particles show larger 
fluctuations, we infer that they may be largely responsible for 
facilitating local rearrangements, in agreement with the find- 
ings of Lynch et al.—. 

X™ 3 *, Af ma x, an d AL max all vary with 0; this dependence is 
shown in Fig. Qj] Both %™ ax and Af max show an increase with 
illustrating that upon approaching the glass transition the dy- 
namic heterogeneity and the associated time scale increases. 
The increasing time scale also suggests that local rearrange- 
ments take longer at higher 0, in agreement with the extended 
plateau at higher in Fig. [3] The characteristic length scale 
AL decreases. This is in excellent agreement with prior work, 
which showed that the displacements for cage rearrangements 
are smaller as the glass transition is approached—. In other 
words, it requires a smaller displacement to be an anomalously 
mobile particle. This can also be seen by comparing Figs. [3] 
and [5] for samples with larger 0, the mean square displace- 
ment has a smaller value even when the non-Gaussian param- 
eter is large, showing that the distribution of displacements is 
overall narrower despite the relatively large fraction of larger- 
than-expected displacements. 

Using the %™ ax data in Fig.[T3]'a) a correlation length scale 
can be estimated by assuming that the correlations %4 mea- 
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Fig. 12 (Color online) (a) is a plot of the susceptibility of large 
particles for various packing fractions, and (b) is a plot of the 
susceptibility of small particles for various packing fractions. 



sures are correlated particles forming compact clusters. Since 
%4 is the average number of correlated particles, then = 
(4/3)71^4, where ^4 is the radius of the cluster of correlated 



particles in units of particle diameters d 



45,64 



The inset in 



Fig. [T3f a) shows the dependence of ^4 on 0. Similarly as 
with the relaxation time, we see a tendency in ^4 to increase 
with 0. The growth in ^4 is about a factor of 4 when the vol- 
ume fraction is increased from a liquid to a dense supercooled 
state. Our values of ^4 are roughly the same as those measured 
in a 2D fluidized granular bed on approaching the jamming 
point—. When compared to i; shown in Fig. [10] the diameter 
of these correlated clusters 2^4 is roughly the same size. 

The time scales Af max are analogous to the a relaxation 
time scales measured in molecular supercooled liquids. In 
many cases the a relaxation time scales are well described us- 
ing either a Vogel-Fulcher-Tammann (VFT) model or Mode- 
Coupling Theory (MCT). 

The first model, VFT, predicts that the time scales should 
obey the form 



AW = Af exp(£/ ( 1 - 0/(|)o) ) , 



(6) 



where A?q, E, and 0o are all fitting parameters. In the model 
Afn is an attempt time to undergo relaxation events over some 
typical length scale. For our experiment, this length scale 
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Fig. 13 (Color online) (a) Plot of the maximum of %4 as a function 
of 0, showing how temporal heterogeneity increases as 
— > g ~ 0.58. The inset in (a) shows the dependence of the 
dynamic heterogeneity length scale ^4 = (x™") 1 '' 3 on 0. (b) Plot of 
the dynamic heterogeneity time scale as a function of 0. (c) Plot of 
the length scale AL as a function of 0. For all panels, the symbols 
are as indicated in the legend of panel (b). 



would be on the order of a particle diameter and the attempt 
time would be the time it takes a particle to diffuse over this 
length scale in the dilute limit. Using the Stokes-Einstein- 
Sutherland formula and a viscosity of 2.18 mPa-s (measured 
for the fluid in absence of colloids) we estimate that at room 
temperature it should take the small particles about 11 sec- 
onds and the large particles 25 seconds to diffuse their own 
diamete r 48 ' 49 . The fitting parameter 0o is the packing fraction 
at which diffusive motion should cease. This should occur at 
random close packing of ~ 0.65 (using the value appropri- 
ate for our binary suspension). However, as pointed out by 
Brambilla et. alJ^, there is a debate as to whether the diver- 
gence predicted by VFT should occur at the jamming point or 
at a slightly different packing fraction. To definitively show if 
this is the case one would need very careful measurements ex- 
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Fig. 14 (Color online) (a) Log-linear and (b) log-log plots of 
temporal dynamical heterogeneity time scale, with fits to eqn ® in 
panel (a) and eqn in panel (b). (c) Log-linear and (d) log-log 
plots of the %4 length scales, with fits to eqn l(6j in panel (c) and 
eqn (|7) in panel (d). 



tremely close to the jamming point which is beyond the scope 
of this paper. The final fitting parameter is E, the fragility, 
which is a material dependent value. The fragility is a mea- 
sure of how sensitive the time scale is to small changes in vol- 
ume fraction. For a molecular system E measures how sen- 
sitive the relaxation time is to small changes in temperature, 
and E ranges between rs 1-100. Materials with low E values 
are termed fragile glass formers and those with large E values 
termed strong glass formers 1 . 

The second model, MCT, predicts a scaling of 



Mnax = Af (l 



(7) 



where Afo, y, and (j) c are the fitting parameters 65 . (j) c in this 
model takes a different meaning with the divergence predicted 
to occur near the glass transition volume fraction, not at ran- 
dom close packing. In light scattering experiments performed 
by Brambilla et. al on 10% poly disperse colloidal samples 
they found (j) c w 0.59, slightly above the glass transition vol- 
ume fraction 14 . Their work also showed that near the diver- 
gence point the dynamics deviate from the predicted form, but 
that in the supercooled regime the MCT equation describes the 
data well. They also found a scaling exponent of y = 2.5 ±0.1. 

Fits to the measured time scales using the two fitting models 
are shown in Fig.[l4fa) and (b), and the corresponding fitting 
parameters are shown in Table Q] In the previous paragraphs 



Table 1 This table displays the fitting parameters found when fitting 
the data to either a VFT scaling or a power law scaling. The 
uncertainties of the fitting parameters are found by adjusting the 
fitting parameters until they no longer provide reasonable fits. 



VFT: 


Af [sec] or Q [d] 


E 


<t>0 


All Af max 


70 ±50 


0.6±0.3 


0.64 ±0.03 


Big A/ max 


200 ±160 


0.4±0.2 


0.6 ±0.03 


Small Af max 


25 ±20 


0.4±0.35 


0.67 ±0.03 


All Ia 


0.3±0.2 


0.4 ±0.2 


0.65 ±0.05 


Big ^4 


0.5±0.2 


0.3±0.2 


0.66 ±0.07 


Small ^ 4 


0.2±0.15 


1.0±0.8 


0.68 ±0.07 




Power law: 


Af [sec] or ^ [d] 


y or 8 




All Af max 


100 ±90 


1.3±0.6 


0.57 ±0.02 


Big Af max 


90 ±70 


1.6±0.8 


0.57 ±0.02 


Small Af max 


20 ±10 


2.6±0.9 


0.61 ±0.01 


All §4 


0.2±0.1 


0.9 ±0.4 


0.59 ±0.03 


Big ^4 


0.5±0.1 


1.4±0.5 


0.62 ±0.04 


Small ^ 4 


0.15 ±0.05 


0.4 ±0.2 


0.57 ±0.02 



reasonable fitting values were given for some of the different 
fitting parameters. The Afo values are significantly larger than 
the dilute concentration diffusive time scales, for both the VFT 
and MCT fits, although the agreement is off by only a factor 
of two for the small particles (25 s for VFT, 20 s for MCT, and 
Id = 1 1 s). For the VFT fit, <|) c is near § rcp as predicted. For 
the MCT fit § c is near the expected glass transition volume 
fraction of m 0.58. The MCT exponent y is smaller than that 
found by Brambilla et. al^, who found y = 2.5, with the ex- 
ception of the small particles for which we find y = 2.6 ± 0.9. 

Our data gives fragilities on the order of 0.5, consistent with 
fragility values from a study of a 2D fluidized granular bed 45 . 
When compared to a molecular system our colloidal system 
would be considered a very fragile glass former. 

In the study on the motion of grains in a 2D fluidized gran- 
ular bed it was shown that the length scales can also be fitted 
well to the models used to fit the time scales where the VFT 
formula becomes ^4 = ^exp(E/(l — t|>/cj> c )) and the MCT 
formula becomes ^4 = c^(l — (j)/(|) c ) 8 , where 8 in work by 
Berthier et. al is predicted to be 2/3 66 67 . The work of Bram- 
billa et. al found that 8 = 2/3 fitted their light scattering data 
very well—. 

The fits to the length scales are shown in Fig. fT4T c) and 
(d), and the fitting values are shown in Table [T] The fitting 
values found for the VFT fits are physically feasible where the 
fragilities and divergence points compare well to the fitting pa- 
rameters previously found for the VFT fits to the time scales. 
The MCT fits are also reasonable, although our scaling expo- 
nents 8 is only consistent with the predicted value of 2/3 due 
to our large error bars. The MCT divergence at <|) e is close to 



10| 



tyg, as expected. 

With the range of volume fractions presented in this paper 
we can not conclusively show which model fits better. Both 
models capture and predict the time and length scales associ- 
ated with dynamic heterogeneity, and the derived fitting pa- 
rameters of both compare well to expected values. 

4 Summary 

We have used confocal microscopy to study three-dimensional 
motion of particles in binary colloidal mixture. The volume 
fraction is varied from 0.4-0.7 and a glass transition, charac- 
terized by aging dynamics, is found at rs 0.58. The dynam- 
ics of large and small particles are qualitatively similar. At 
volume fractions approaching the glass transition, both show 
an increase in motion at the same characteristic cage break- 
ing time scale. This time scale also corresponds with the time 
over which the displacement distribution functions are broad- 
est (most non-Gaussian). Particle motion is facilitated by the 
presence of small neighbors, and inhibited by large neighbors, 
consistent with the idea that small particles serve as lubricants. 
We have investigated vector and scalar correlation functions 
and extracted specific length scales associated with the spatial 
decay in correlation of the displacements. This length slightly 
increases with volume fraction, although it does not appear 
to diverge as the glass transition is approached. The tempo- 
ral correlations also give rise to length scales and time scales 
which grow as the glass transition is approached, although the 
form of this growth is ambiguous with respect to power-law 
or exponential growth. 

The presence of particle tracking noise makes certain mea- 
surements more difficult, in particular, the pair correlation 
functions (Pig. [2j and the mean square displacement (Fig . [3} . 
The primary conclusions of our work, however, focus on 
the particles that move large distances, and these measure- 
ments have a "signal" (the distance moved) larger than the 
"noise" (the instantaneous positional uncertainty). Our mea- 
surements of the non-Gaussian parameter (Pig. [5] spatial cor- 
relation functions (Figs. [8l9T >. correlation lengths (Fig.[l0] and 
dynamic susceptibility measurements (Figs. QT| - [T4T i are ro- 
bust to the noise. Likewise, the identification of nearest neigh- 
bors is fairly robust to even moderate fluctuations in pair-wise 
particle separations, and so we have confidence in our data 
showing that having fewer large neighbors enhances a parti- 
cle's mobility (Fig. |6j. 
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